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Abstract 

This paper describes an integrated aerodynam- 
ic/dynamic/structural (IADS) optimization procedure for 
helicopter rotor blades. The procedure combines perfor- 
mance, dynamics, and structural analyses with a general 
purpose optimizer using multilevel decomposition tech- 
niques. At the upper level, the blade structure and re- 
sponse are represented in terms of global quantities 
(stiffnesses, mass, and average strains). At the lower level, 
the blade structure and response are represented in terms of 
local quantities (detailed dimensions and stresses). 

The upper level objective function is a linear combina- 
tion of performance and dynamic measures. Upper level 
design variables include pretwist, point of taper initiation, 
taper ratio, root chord, blade stiffnesses, tuning masses, 
and tuning mass locations. Upper level constraints consist 
of limits on power required in hover, forward flight, and 
maneuver; airfoil drag; minimum tip chord; trim; blade 
natural frequencies; autorotational inertia; blade weight; 
and average strains. 

The lower level sizes the internal blade structure at sev- 
eral radial locations along the blade. The lower level op- 
timization assures that a structure can be sized to provide 
the stiffnesses required by the upper level and assures the 
structural integrity of the blade. The lower level design 
variables are the box beam wall thicknesses and several 
lumped areas that are analogous to longitudinal stringers in 
a wing box cross section. The lower level objective func- 
tion is a measure of the difference between the upper level 
stiffnesses and the stiffnesses computed from the wall 
thicknesses and lumped areas. Lower level constraints are 
on the Von Mises stress at the box comers for multiple 


load cases generated by several flight conditions, limits on 
wall thicknesses for thin wall theory, and other dimen- 
sional considerations. 

The IADS procedure provides an optimization tech- 
nique that is compatible with industrial design practices in 
which the aerodynamic and dynamic design is performed 
at a global level and the structural design is carried out at a 
detailed level with considerable dialogue and compromise 
among the aerodynamic, dynamic, and structural groups. 
The IADS procedure is demonstrated for several cases. 

Notation 

A area (ft 2 ) 

AI autorotational inertia, ^w r? (lbm-ft 2 ) 

j=i 

a i ith lumped area (ft 2 ) 

b box width (ft) 

c d airfoil section drag coefficient 

Cd^j maximum allowable section drag coefficient 
c Y largest section drag coefficient at azimuthal angle w 

°max 

CF centrifugal force Ob) 

C D rotor coefficient of drag 

C L rotor coefficient of lift 

c f root chord (ft) 

c ( tip chord (ft) 

DV pth upper level design variable 

E Young’s modulus of elasticity (lb/ft ) 

EA extensional stiffness (lb) 

EI XX chordwise bending stiffness (lb-ft 2 ) 

EI^ flapwise bending stiffness (lb-ft 2 ) 
f b- ith bending frequency (per rev) 



f k kth frequency (per rev) 

f u lower bound on kth frequency (per rev) 
r ku upper bound on kth frequency (per rev) 

f t j itfa torsional frequency (per rev) 

F lower level objective function 

g Cj ith lower level constraint function 

g m „ maximum lower level constraint function, 
max{ g c . } 

ith upper level constraint function 

G torsional modulus erf elasticity (lb/ft 2 ) 

GJ torsional stiffness (lb- ft 2 ) 

h box height (ft) 

I xx chordwise moment of inertia (ft 4 ) 

1^ flapwise moment of inertia (ft 4 ) 

ITER number of trim iterations 
J polar moment of inertia (ft 4 ) 

kj ith weighting factor in objective function 

KS Kresselmeir-Steinhauser function 

If factor of safety 

mj ith segment tuning mass (slug/ft) 

N number of blades 

n integer 

nc number of constraints components in lower level 
nP frequency at n times rotational speed of blade 
NDV number of upper level design variables 

OBJ upper level objective function 

P main rotor power (hp) 

R blade radius from center of rotation (ft) 

r distance along blade from center of rotation (ft) 

Tj distance from center of rotation to center of jth 

segment (ft) 

$Nff N per rev vertical rotating hub shear in forward 
flight (lbf) 

J> N ref reference N per rev rotating vertical hub shear in 

forward flight (lbf) 
t k kth wall thickness (ft) 

v A ith lower level design variable 

V(o,x) Von Mises stress (lb/ft 2 ) 

W total blade weight (lbm) 

Wj total weight of jth structural segment (lbm) 
y- location of ith tuning mass 

y^ point of taper initiation 

Af increment used in frequency window (per rev) 

e coordination parameter 

ty average strain 

X Lagrange multiplier 

p pull down factor 

o bending stress (lb/ft 2 ) 

x shear stress (lb/ft 2 ) 

¥ azimuth angle, zero over tail (deg) 


0 tw maximum pretwist (deg) 
Subscripts 

a available or allowable 

ff forward flight 

h hover 

m maneuver 

max maximum 

min minimum 

ref reference 

Superscripts 
L lower level 

U upper level 

nondimensional 


Introduction 

Over the last decade optimization techniques have been 
studied for application to the rotor blade design process. In 
Ref. 1 Miura presents a survey chi the application of nu- 
merical optimization methods to helicopter design prob- 
lems including rotor blade design. Most optimization pro- 
cedures have dealt with a single discipline such as aerody- 
namics (Refs. 2-4), structures (Ref. 5), or dynamics (Refs. 
2, 6-9). However, the rotor blade design process is multi- 
disciplinary involving couplings and interactions between 
several disciplines such as aerodynamics, dynamics, struc- 
tures, and acoustics. These couplings and interactions can 
be exploited by the optimization procedure if all the disci- 
plines are accounted for simultaneously rather than sequen- 
tially. For instance, in a review (Ref. 10) on the impact of 
structural optimization on vibration reduction, Friedmann 
emphasizes the need to include the multidisciplinary cou- 
plings between aerodynamics, dynamics, and structures 
even when optimizing only for minimum vibration. 

Techniques and strategies for merging disciplines to ob- 
tain integrated rotorcraft optimization procedures are de- 
veloping. In Refs. 1 1 and 12, a plan is described for inte- 
grating the disciplines of aerodynamics, dynamics, struc- 
tures, and acoustics. As part of that plan, aerodynamics 
and dynamics have been incorporated systematically into 
performance (Refs. 3 and 4) and airload/dynamic (Ref. 13) 
optimization procedures resulting in an integrated aerody- 
namic/dynamic optimization procedure (Ref. 14). Ref. 15 
summarizes recent accomplishments based on that plan. 

Other multidisciplinary rotor blade optimization work is 
described in Refs. 16-19. Refs. 16 and 17 describe the 
formulation of a multidisciplinary approach to rotor blade 
design for improved performance and reduced fuselage vi- 
brations. Ref. 18 describes a staged optimization proce- 
dure for a rotor for combined aerodynamics, dynamics, and 
structures. Ref. 19 describes a multidisciplinary optimiza- 
tion procedure to design high speed prop rotors. 
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What is lacking in previous multidisciplinary rotor 
blade optimization procedures is an efficient method to in- 
tegrate structures or structural properties. Usually struc- 
tures or structural properties are included in one of two 
ways - either as local design variables (indirectly affecting 
the response of the blade) or global design variables 
(directly affecting the response of the blade). When local 
design variables are used, the detail dimensions of a struc- 
tural member at one or more radial locations along the 
blade are used to generate structural properties. When 
global design variables are used, structural properties are 
the design variables. Both type of design variables have 
limitations. Using local design variables (e.g., Refs. 
6,7,18-19), such as wall thicknesses of the structural mem- 
ber, can lead to a large number of design variables which 
can be computationally expensive. Also, this choice of de- 
sign variables is at odds with traditional design practice 
where chord, stiffness, and mass distributions along the 
blade are determined and then a structure is designed 
which matches these distributions. Using global design 
variables (e.g., Refs. 2,9,13,14,16-17), such as stiffness and 
mass properties, in optimization also has disadvantages. 
When flapwise bending stiffness, chordwise bending stiff- 
ness, torsional stiffness, and extensional stiffness distribu- 
tions are used as design variables, they are treated as inde- 
pendent quantities. In reality, these stiffnesses are not in- 
dependent, and there is no guarantee that a set of wall 
thicknesses can be found which will simultaneously give 
these stiffnesses. 

This paper presents the methodology for incorporating 
aerodynamics, dynamics, and structures in an integrated 
optimization procedure using both local and global design 
variables. Multilevel decomposition techniques based on 
Ref. 20 are used to add structural design variables and 
constraints to an existing aerodynamic/dynamic optimiza- 
tion procedure (Ref. 14). The product is an integrated 
aerodynamic/dynamic/structural optimization (IADS) pro- 
cedure. The multilevel decomposition formulation used in 
this paper was presented first in Ref. 15. Another prelimi- 
nary study of multilevel decomposition techniques applied 
to rotor blade design is described in Ref. 21. 

The multilevel decomposition approach has been suc- 
cessfully applied to multidisciplinary problems (e.g., Refs. 
22-24). As originally proposed in Ref. 25, the coordination 
procedure consisted of an optimum sensitivity analysis 
(Ref. 26) and a set of equality constraints which relate the 
detailed (local) design variables of one subsystem to the 
global design variables on the level above. However, as 
pointed out in Ref. 27, these equality constraints have 
caused difficulties in implementing multilevel decomposi- 
tion procedures. The IADS procedure is based on the mul- 
tilevel decomposition approach of Ref. 20 which elimi- 
nates the equality constraints in the coordination procedure 
allowing the use of the less computationally costly opti- 
mum sensitivity derivative found in Ref. 28. However in 
the IADS procedure, the set of lower level constraints is 


replaced by an envelope function known as the 
Kresselmeir-Steinhauser function (KS-function, Ref. 29) 
which further reduces computational cost 

First, the general multilevel decomposition strategy 
with two levels will be discussed (note: systems with more 
levels are discussed in Refs, 20, 22, and 25). Next, the 
general strategy will be related to rotor blade design. 
Then, the IADS development including flowcharts of the 
upper and lower levels and the optimization procedure will 
be explained. Results will be presented for several cases 
which demonstrate the strengths of the IADS procedure. 


Multilevel Decomposition Optimization Strategy 

With a multilevel decomposition approach (Refs. 20, 
22, and 25), a large complex optimization problem is bro- 
ken into a hierarchy of smaller optimization subproblems. 
This hierarchy can be thought of as levels of increasing 
detail. At the upper level, the subproblem is formulated in 
terms of global quantities which describe the overall be- 
havior of the entire system. On the lower level, the sub- 
problems are stated in terms of local quantities and local 
constraints which have only a small impact on the entire 
system. Each subproblem uses local design variables to 
reduce the violation of constraints which are unique to that 
subproblem. The coupling between the upper level sub- 
problem and the lower level subproblems is preserved 
through a coordination procedure such as that described in 
Refs. 20 or 25. This coupling represents a dialogue be- 
tween the levels that upon convergence establishes compat- 
ibility between the two levels. 

Fig. 1 illustrates a generic two-level optimization pro- 
cedure. Note that the analysis proceeds from the upper 
level to the lower level while the optimization proceeds 
from the lower level to the upper level. First, the upper 
level analysis initializes all the global quantities and re- 
sponses and then provides information to each lower level 
subproblem. Individual lower level optimizations are per- 
formed which reduce local constraint violations as much as 
possible and which provides information to the coordina- 
tion procedure. Next the upper level optimization occurs. 
The preceding describes one cycle. The entire process is 
repeated for several cycles. Convergence occurs when all 
the constraints (both upper level and lower level) are satis- 
fied and the upper level objective function is minimized. 

The rotor blade optimization problem can be decom- 
posed into one subproblem affecting the global response of 
the blade and three subproblems affecting portions of the 
blade. Quantities such as power required, blade trim, au- 
torotational inertia, natural frequencies, total blade weight, 
and average strain describe the global response of the 
blade. The entire blade must be analyzed to obtain these 
response quantities. Quantities such as stresses are detailed 
response quantities since only a portion of the blade must 
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be analyzed to obtain these response quantities. Therefore, 
a two-level decomposed rotor blade optimization problem 
can be defined as shown in Fig. 2. The upper level 
optimizes the blade by changing global quantities such as 
blade planform, twist, and distributions of mass and 
stiffness. The upper level chord, mass, and stiffness 
distributions are treated as independent quantities. The 
reconciliation between these distributions is done on the 
lower level. The lower level consists of several indepen- 
dent subproblems at stations along the blade radius, which 
optimize detailed cross-sectional dimensions to satisfy 
stress constraints and reconcile the upper level independent 
mass, chord, and stiffness distributions with the lower level 
calculated stiffness distributions. This reconciliation is 
improved further by a set of upper level coordination con- 
straints (see Appendix A). First the upper level analysis 
and optimization will be described, then the lower level 
analysis and optimization, and then the overall IADS 
system. 


Upper Level Analysis and Optimization 

The purpose of the upper level analysis is to evaluate 
the overall rotor blade design on the basis of performance, 
dynamic, and global structural measures. (For a descrip- 
tion of the rotor blade design philosophy see Refs. 3-4, 11- 
12, and 14-15.) The upper level analysis is similar to the 
integrated aerodynamic/dynamic analysis reported in Ref. 
14 with the addition of extensional stiffness design vari- 
ables, strain constraints, and coordination constraints. As 
shown in Fig. 3, the blade is analyzed for three flight con- 
ditions: hover, forward flight, and maneuver. The 
Langley-developed hover analysis program HOVT (a blade 
element momentum analysis based on Ref. 30) is used to 
predict power required in hover. The comprehensive heli- 
copter analysis program CAMRAD/JA (Ref. 31) is used to 
predict rotor performance (e.g., trim, airfoil drag, power 
required), loads, and frequencies for forward flight and 
maneuver. The maneuver flight condition simulates a co- 
ordinated turn in terms of a increased load on the forward 
flight lift requirement 

Upper Level Design Variables - The upper level design 
variables are the blade planform, stiffnesses, and tuning 
masses (see Fig. 4). The blade planform is defined by the 
point of taper initiation y^, root chord c r , taper ratio c^, 

and maximum pretwist 0 tw . The blade is rectangular from 
the root to y^ and then tapers linearly to the tip. The 
pretwist varies linearly from the center of rotation to the 
tip. Global design variables include the blade chordwise, 
flapwise, torsional, and extensional stiffnesses (denoted by 
EI XX , EI^ GJ, and EA, respectively) at three radial loca- 
tions: blade root, point of taper initiation, and blade tip. 
The stiffnesses are assumed to vary linearly between these 
points and are treated as independent quantities. The re- 
maining design variables are three tuning masses (denoted 


by mj, mj, m 3 ) and their locations (denoted by y y 2 , and 
y 3 ). The total blade mass consists of the structural mass 

(which is assumed constant) plus the sum of the tuning 
masses. There is no attempt to reconcile the change in 
weight with the change in design variables since the pre- 
sent work is based on extending the procedure of Ref. 14 to 
include structures. However, this reconciliation is possible 
(see Ref. 15). It is assumed that the center of gravity and 
aerodynamic offsets are coincident with the blade elastic 
axis. The number of blades, rotor radius, rotational veloc- 
ity, airfoils, and airfoil distribution are preselected and 
fixed. 

Upper Level Objective function - The objective function 
to be minimized is a combination of performance and dy- 
namics measures and is formulated as follows 


OBJrrkj 




r m -+*4 

m ref 


lUSL 

S N ref 


( 1 ) 


where p h’ p ff and P m are the powers required in hover, 

forward flight, and maneuver, respectively. N is the num- 
ber of blades and is the N per rev rotating vertical 

hub shear in forward flight. The terms kj, k 2 , k 3 , and k 4 
are weighting factors chosen by the user. P href , Pff rcf > 
Pm rc f » and SN ref are reference values used to normalize 

and nondimensionalize the objective function components. 
The usefulness of this objective function was demonstrated 
in Ref. 14. 

Upper Level Constraints - The rotor blade design pro- 
cess is defined in terms of aerodynamic performance, dy- 
namics, and global structural requirements. Satisfactory 
aerodynamic performance is defined by the following four 
requirements. First, the power required for any flight 
condition must be less than the available power. Second, 
airfoil section drag along the blade radius on the advancing 
and retreating side of the rotor disk in both forward flight 
and maneuver must be less than a maximum allowable 
value. Third, the rotor must trim at each flight condition. 
The rotor is trimmed to a constant lift in forward flight and 
a (different) constant lift in maneuver which ensures that 
the rotor has no loss in lift capability or maneuverability 
even if solidity decreases from the initial to the final de- 
sign. Incorporation of a maneuver flight condition is used 
in place of a constraint on solidity, since low speed ma- 
neuver determines rotor solidity (Ref. 32). Fourth, the 
blade tip chord must be larger than a prescribed minimum 
value. Satisfactory dynamics is defined in terms of limits 
on vibrational frequencies. The blade is designed so that 
the natural frequencies (both bending and torsional) do not 
coincide with integer multiples of the rotor speed. Also, 
the blade must have sufficient autorotational inertia as a 
safety measure needed in case of engine failure. In addi- 
tion to satisfying these design requirements, the blade 
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weight must not exceed some upper limit. Satisfactory 
structural requirements are defined in terms of limits on the 
average axial strains for forward flight and the maneuver 
flight conditions. 

This section of the paper discusses the performance, 
dynamic, and structural constraints. The coordination con- 
straints are discussed later in the paper. The performance 
and dynamic constraints are the same as those used in Ref. 
14. By convention, the ith constraint gj is satisfied if it is 
less than or equal to zero. 

Performance Constraints - The performance con- 
straints are on power required, trim, airfoil section drag, 
and blade tip chord. The requirement that the power re- 
quired be less than the power available is given by 

p i 

g i = — - -1^0 for each flight condition (2) 

*a 

where Pj is the power required for the ith flight condition 
and P a is the power available. 

The requirement on the airfoil section drag translates 
into a constraint that each airfoil section distributed along 
the rotor blade operate at a section drag coefficient c^ less 

than a specified allowable value c<j al| (see Ref. 14). This 

leads to 24 constraints per flight condition since the blade 
is analyzed in 15 degree azimuthal increments around the 
rotor disk. At a given azimuthal angle Y the constraint is 
formulated as 

c* 

g . = max_ _ l <; 0 «j/ =15 3Q 45 m (3) 

Cd »ll 

Uf 

where c<j is the allowable drag coefficient and cj is 

ail a max 

the largest drag coefficient at any radial station (note: the 
drag coefficients in the reverse flow region occurring on 
the retreating side of the rotor disc are ignored). In the pre- 
sent work the same value for Cd a|| is used on the advanc- 
ing and retreating side of the rotor disk. This simplifying 
assumption could easily be lifted. 

The trim requirement is difficult to translate into a 
mathematical constraint. The trim constraints in forward 
flight and maneuver are implemented using the method de- 
veloped in Ref. 3 which expresses the constraint in terms 
of the number of trim iterations (ITER), the maximum 
number of trim iterations allowed (ITER max ), and the pth 

nondimensional design variable ( DVp). The heuristic trim 
constraint is given by 


NDV 

gi = (ITER - ITER max + 1)( X DV p ) £ 0 (4) 

P=1 

where NDV is the number of design variables. In devel- 
opment of this equation in Ref. 3, it was found that the 
addition of the summation term improved convergence be- 
cause it allowed calculation of the change in the trim con- 
straint with respect to change in a single design variable. 

The final performance requirement is a constraint used 
to ensure that the blade tip chord does not become too 
small 

gi =1 — ^0 (5) 

Cf • 
l min 

where c f is the tip chord and c t is the minimum tip 

l l min r 

chord allowed. This is a practical constraint used to assure 
validity of the airfoil tables and address manufacturing 
considerations. 

Dynamic Constraints - The dynamic constraints are 
on frequencies, total blade weight, and autorotational iner- 
tia. The constraint on the kth frequency f^ (either a bend- 
ing or a torsional frequency) is formulated such that the 
frequency is separated from integer multiples of the rotor 
speed by an amount Af 


o 

VI 

7 

n 

& 

(upper bound) 

(6a) 

*ku 



and 



gi =l--p L <0 

f kl 

(lower bound) 

(6b) 

where f ku has a value lhat is Af below n+1 per rev and f k( 


has a value that is Af above n per rev for the applicable n. 
For example, suppose Af is 0.1 per rev and f 4 is 5.6 per 

rev, then nP would be 5 per rev and (n+l)P would be 6 per 
rev. Thus f 4u and f 41 would be 5.9 per rev and 5.1 per rev, 

respectively. Formulating the constraints in this manner 
allows the frequencies to change from one optimization 
cycle to the next cycle provided the frequencies avoid ap- 
proaching integer multiples of the rotor speed. This for- 
mulation is different from the approaches used in Refs. 13 
and 16-17 where the frequencies are kept within prescribed 
windows based on the reference blade frequencies. In this 
work, constraints are placed on frequencies in both forward 
flight and maneuver since blade collective pitch and the 
amount of modal coupling may be different for the two 
flight conditions and therefore the frequencies could be 
different. 

The constraint that the blade weight be less than some 
maximum value is formulated as follows 


5 



w 


-ISO 


( 12 ) 


8i = 


W„ 


(7) 


NDV 
g = go + Z 


iL 


M ®V { 


ADVj 


where W is the total blade weight and W max is the maxi- 
mum allowable weight. The total blade weight is the 
structural mass distribution (which is constant) plus the 
sum of the tuning masses. 

Finally, the blade must have enough autorotadonal iner- 
tia AI for safe autorotation in case of engine failure. The 
constraint is formulated so that the autorotational inertia of 

the blade is greater than some minimum value AI . 

min 

gi = l--7^-*0 (8) 

A^min 


Structural Constraints - The structural constraints are 
on the average axial strains. The structural constraints 
which are evaluated at the same radial locations used to 
define the design variables (Fig. 4) are imposed on the av- 
erage axial strains e as follows 


and 


E Y 

g i =- 2 --iso 

Ey 

£ a 


(9a) 

(9b) 


where e a is (he magnitude of the allowable strain and 


e 


y - 


IjC F 
EA 


(10) 


where CF is the centrifugal force, EA is the extensional 
stiffness, and If is a safety factor on the loads. The strain 
constraints are calculated using loads from both the for- 
ward flight and the maneuver flight conditions. 

Upper Level Optimization - The upper level optimiza- 
tion consists of the general purpose optimization program 
CONMIN (Ref. 33) and an approximate analysis used to 
reduce the number of HOVT and CAMRAD/JA analyses 
during the iteration process. The approximate analysis is 
used to extrapolate the upper level objective function and 
upper level constraints with linear Taylor Series expan- 
sions using derivatives of the objective function and con- 
straints with respect to the design variables 

NDV 

OBJ = OBJ 0 + 2 ADVj (11) 

i=l ” DV i o 


The assumption of linearity is valid over a suitably small 
change in the design variable values and will not introduce 
a large error into the analysis provided the changes ADV* 
are small. Errors which may be introduced by use of the 
approximate analysis are controlled by imposing “move 
limits*' on each design variable during the iteration process. 
A move limit which is specified as a fractional change of 
each design variable value is imposed as an upper and 
lower design variable bound. At the present time the move 
limits are manually adjusted. 


Lower Level Analysis and Optimization 

This section of the paper describes the lower level anal- 
ysis and lower level optimization procedure. The puipose 
of each lower level optimization is to assess whether a 
structure at the given radial location can be sized to pro- 
vide the stiffnesses required by the upper level optimiza- 
tion and have the strength to withstand loads calculated by 
the upper level analysis. The lower level optimizations can 
be done in parallel since they are independent. 

For simplicity, since closed-form equations can be de- 
rived (see Appendix B), the structural member (Fig. 5) is 
assumed to be a thin-walled isotropic box. The box cross 
section is symmetric about the horizontal axis with wall 
thicknesses ^ and lumped areas aj which are analogous to 
longitudinal stringers in a wing box cross section. The 
outer dimensions b (the box beam width) and h (the box 
beam height) are functions of the upper level design vari- 
ables since b and h depend on the local chord and the local 
airfoil thickness. The values of b and h are determined by 
placing a box of maximum area within the airfoil cross 
section using the method of Ref. 34. 

Lower Level Design Variables - The design variables are 
the three wall thicknesses (tj, tj, and t 4 ) and the three 

lumped areas (a^ ^ and a 3 ). The lumped areas are used 
to give the lower level more flexibility in matching the 
upper level stiffnesses. For the present implementation, 
the lumped areas are assumed to be square areas. 

Lower Level Objective Function - The objective func- 
tion is a measure of the difference between the stiffnesses 
required on the upper level and those determined from the 
lower design variables 
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where a starred quantity ( )* denotes an upper level design 
variable. The lower level cross sectional properties I xx , 
Ijx and J are computed (see Appendix B), E is Young's 
modulus of elasticity, and G is the torsional modulus of 
elasticity. 

Lower Level Constraints - The constraints are enforced 
on the extensional stiffness, stresses, and the physical di- 
mensions of the wall thicknesses and lumped areas. The 
extensional stiffness constraint which requires the lower 
level calculated extensional stiffness EA (see Appendix B) 
to be equal or greater than the upper level extensional stiff- 
ness (EA)* (an upper level design variable) is given by 


g Ci =1“ 


EA 

(EA)* 


<0 


(14) 


at the given cross section. It is noted that the extensional 
stiffnesses appear in a constraint rather than in the objec- 
tive function (Eq. 13) where the other stiffnesses appear. 
This is done for the following reason. The role of EA in the 
upper level is limited to satisfying the strain constraints 
(Eq. 9). The lower level is responsible only for assuring 
that the value of EA is at least as large as the value needed 
in the upper level - close matching of EA to (EA)* is not 
required. 

The stress constraint which is evaluated at the comer of 
the box cross section shown in Fig. 5 has the following 
form 


8e =^(0,1) 
&c i rr 


(15) 


where a is the bending stress, x is the shear stress, and 
V(o,x) is the Von Mises stress measure (see Appendix B). 
Two stress constraints are used - in one x is based on the 
vertical wall thickness and in the other x is based on the 
horizontal wall thickness. 

A set of constraints is imposed on the lower level wall 
thicknesses to assures that the section remains a thin- 
walled section and that the expression for J remains valid 


(see Appendix B). These constraints are 


ti 

g c = — 1<0 

c ‘ 0.1b 

j = 2 and 4 

(16) 

t; 

g c =~ J 1^0 

6C ' O.lh 

j = 1 and 3 

(17) 


where b and h are the width and height of the box cross 
section, respectively. 


A set of constraints is imposed on the lumped areas and 
wall thicknesses which require that the dimensions are 
physically possible (i.e. the lumped areas can Fit inside of 
the box cross section. These constraints are 

= _[b-^_^_^] so (lg) 

8c, (19) 


g Cj = -[h - 1, - 1 3 - 2^] £ 0 (20) 

g Cj = -[h - tx - 1 3 - 2yfe] < 0 (21) 

g Cj = -[h - ti - 1 3 - 2yfe] < 0 (22) 

In addition a set of constraints representing upper and 
lower bounds on the design variables is used. For the kth 
design variable, the lower bound is given by 

g C j = Vkj - v k :£ 0 (23) 

and the upper bound by 

gq = v k - v k u - 0 (24) 


where v ^ and vj^ are the lower and upper design vari- 
able bound, respectively. 

For convenience, the set of lower level constraints de- 
fined by Eqs. 14-24 is replaced by a single cumulative con- 
straint, an envelope function known as the KS function 
(Ref. 29), which approximates the active constraint bound- 
ary 


KS — gmax 

P 


P(8c; 6 max) 

Xe J 

j=l 


<0 


(25) 


where g max is the maximum constraint component from 

Eqs. 14-24, nc is the number of lower level constraint com- 
ponents and p is defined by the user. Initially p is small 
and then increases until a maximum value p max is reached. 
For large values of p, the value of KS approaches g max . 

The KS function is a single measure of the degree of con- 
straint satisfaction or violation and is positive (violated) if 
at least one of the constraints g c is violated. The KS 

function is a single-valued function which is continuous 
and differentiable. This property becomes important when 
implementing the upper and lower levels as described in 
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the section on the overall organization of the IADS proce- 
dure. 

Lower Level Optimization Procedure - The flowchart 
for each lower level optimization procedure is shown in 
Fig. 6. Loads, local chord, box beam width, box beam 
height, and upper level stiffnesses are passed down from 
the upper level analysis. The lower level design variables 
(Fig. 5) are used to calculate lower level stiffnesses. Von 
Mises stresses are calculated using the loads from the for- 
ward flight and maneuver analyses. The lower level ob- 
jective function (Eq. 13) and cumulative constraint (Eq. 
25) are evaluated. The lower level optimizations are per- 
formed using the general purpose optimization program 
CONMIN. Exact analyses are used to evaluate the objec- 
tive function, the constraint, and any gradients computed 
by CONMIN. The optimization process is converged 
when the objective function is minimized and the cumula- 
tive constraint is satisfied. After convergence the process 
returns to the upper level. 


Coordination Between Upper and Lower Levels 

The coordination between upper and lower levels is im- 
plemented by upper level constraints. These constraints 
are imposed to encourage changes in the upper level design 
variables which promote consistency between the upper 
and lower level stiffnesses. Specifically, these constraints 
(one for each lower level optimization) have the form 

g = F u -(l + e)F^- 50 (26) 

where F^ is the most recent value of the lower level ob- 
jective function (i.e., optimum value of Eq. 13), F u is an 
estimate of the change in which would be caused by a 

change in the upper level design variable values, and e is a 
specified tolerance denoted the coordination parameter (see 
Appendix A). The importance of this parameter will be 
discussed later. 

Eq. 26 is the general form of the coordination constraint 
as formulated in Ref. 15. As shown in Appendix A, the 
coordination constraint can be approximated in terms of 
the lower level total optimum sensitivity derivative which 
expresses how the optimum lower level objective function 
and lower level active constraint will change with a change 
in upper level design variable. 


Overall Organization of IADS Procedure 

The conceptual IADS procedure is shown in Fig. 2. It 
consists of an upper level analysis (Fig. 3), three lower 
level optimizations (Fig. 6), and a coordination task. The 
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actual IADS procedure is more complicated and requires 
an upper level sensitivity analysis and three lower level op- 
timum sensitivity analyses in addition. 

The flowchart for the IADS procedure is shown in Fig. 
7. First, the upper level analysis is executed for the current 
set of design variables providing all of the information 
needed to calculate the upper level objective function and 
constraints with the exception of the coordination con- 
straints. The upper level analysis also provides the loads, 
local chord, box beam width, box beam height, and stiff- 
nesses (to be matched) to the lower level analysis. F ach 
lower level optimization is performed to obtain a set of 
lower level design variables which match the current upper 
level bending and torsional stiffnesses as close as possible. 

Next, an upper level sensitivity analysis is performed 
consisting of forward finite difference derivatives (or gra- 
dients) of the upper level analysis. These derivatives of the 
upper level objective function and upper level constraints 
are required to approximate the upper level objective func- 
tion and upper level constraints during the upper level op- 
timization. In addition, the loads and local chords corre- 
sponding to the changes in the upper design variables are 
saved. These quantities are used in the three lower level 
optimum sensitivity analyses to approximate the coordina- 
tion constraint (Eq. 26). Appendix A describes how the 
coordination constraint is expressed in terms of the total 
optimum sensitivity derivative involving changes in the op- 
timum lower level objective function with respect to 
changes in the upper level design variables and changes in 
the active lower level constraint with respect to changes in 
the upper level design variables. 

Finally, the upper level optimization consisting of 
CONMIN and approximate analysis occurs. This describes 
one cycle of the IADS procedure. The process is repeated 
for additional cycles until convergence is achieved. A very 
strict convergence criterion is used for demonstration pur- 
poses. The overall procedure is converged when the 
change in the upper level objective function is less than 

0.5x10 over three consecutive cycles and all the con- 
straints (both upper and lower level) are satisfied. A step 
size of 0.001 is used to compute the finite difference 
derivatives. 


Demonstration of the IADS Procedure 

This section of the paper describes the analytical blade 
model, the mission definition, the optimization problem, 
and optimization results used to demonstrate the IADS 
procedure. Results are presented for two studies - (1) the 
effect of initial design and (2) the effect of the coordination 
parameter e. 

Analytical Blade Model - The analytical blade model 
used to demonstrate the IADS procedure represents a wind 



tunnel model of a rotor blade for a four-bladed helicopter. 
The blade has radius is 4.68 ft. Three sets of advanced air- 
foils are used along the blade - the RC(4)-10 airfoil (Ref. 

35) from the root to 85 percent radius, the RC(3)-10 (Ref. 

36) airfoil from 85 to 95 percent radius, and the RC(3)-08 
(Ref. 36) airfoil from 95 percent radius to the tip. Tables 
of experimental two-dimensional airfoil data for these three 
airfoil types are used in both HOVT and CAMRAD/JA. 
The analytical model of the blade uses 19 aerodynamic 
segments for HOVT, 50 structural segments and 18 aero- 
dynamic segments for CAMRAD/JA. HOVT is used to 
predict the power required in hover using nonuniform in- 
flow (no wake is included) by trimming to a constant lift 
Cl* CAMRAD/JA is used to predict rotor performance, 
loads, and frequencies using uniform inflow with empirical 
inflow correction factors for the forward flight and maneu- 
ver flight conditions. Uniform inflow is used to save on 
computational costs (note: even though approximate anal- 
ysis is used in the upper level optimization 46 
CAMRAD/JA analyses are required per optimization cy- 
cle). In CAMRAD/JA an isolated rotor analysis is used 
which trims the rotor to constant lift C L and drag C D and 
zero flapping angle relative to the shaft using collective, 
lateral cyclic, and longitudinal cyclic pitch. From the 
modal analyses in CAMRAD/JA using ten bending modes 
and five torsional modes, it is found that only the first six 
bending frequencies are below 10 per rev and need to be 
constrained for a four-bladed rotor. Since f^ corresponds 

to a rigid body mode and f *> 2 is the 1 per rev, the first two 

frequencies are not constrained. Constraints are placed on 
the first four bending frequencies (fb 3 » fb 4 » a n d fb 6 

flapping-dominated and fb 5 lead-lag dominated) and the 
first two torsional frequencies ( f tj representing the rigid 

body torsional mode due to the control system stiffness and 
f t2 representing the first elastic torsional mode). 

Mission Definition - The flight conditions are a constant 
lift of 1-g (331 pounds, C L =0.0081), propulsive force of 32 

pounds (C D =-0.000811), and an advance ratio of 0.35 for 

the forward flight condition and a constant lift of 401 
pounds (C l = 0.00985), a propulsive force of 23 pounds 

(C D =-0 .000596), and an advance ratio of 0.3 for the ma- 
neuver flight condition. The maneuver flight condition is 
for a load factor of 1 .22. These flight conditions and load 
factor are similar to those used in Ref. 37. 

Optimization problem - The objective function is a com- 
bination of the power required in hover, forward flight, and 
maneuver and the 4 per rev rotating vertical hub shear in 
forward flight. The objective function is chosen to be one 
dominated by performance with little emphasis on dynam- 
ics. Of the three powers it is assumed that it is most im- 
portant to reduce the power required in hover - it will have 
twice the weight as the other two powers. Several values 


were tried for the weighting factor on the hub shear term. 
It was found that to obtain the proper balance between per- 
formance and dynamics, k 4 has to be between one and two 

orders of magnitude less than k ^ . Thus, for this case, the 
weighting factors are chosen to be kj =10.0, k 2 =k 3 =5.0, 
and k 4 = 0.5. 

OBJ = 10-^- + 5^ff_ + 5^EL_ + o.5-|^- (30) 

p h ref P ff re f P m ref S 4 ref 

W here P href , P ffref . P mrcf . and S^are 15 hp. 13 

hp, 12 hp, and 2 lbf, respectively. The reference values are 
chosen to be representative of the powers required and hub 
shear for all the initial blade designs used in this work. 

The upper and lower bounds for the design variables are 
given in Table 1. On the upper level twenty-two design 
variables and 95 constraints are used. On the lower level 
six design variables and one cumulative constraint (the KS- 
function with 24 components) are used at each of the three 
span wise locations (i.e., the root, the point of taper initia- 
tion, and the tip). 

Parameters and flight conditions are summarized in 
Table 2. Since the blade is made of aluminum, E has a 
value of 15.26xl0 8 lb/ft 2 , the allowable strain e a has a 
value of 0.05 ft/ft, and the allowable stress a a is 8.352xl0 6 
lb/ft 2 . The values for minimum tip chord C f . , power 

l min 

available P a , minimum autorotational inertia, and maxi- 
mum allowable drag coefficient are 0.083 ft, 20 hp, 

23.69 lbm-ft 2 and 0.12, respectively. Frequencies must be 
at least 0.1 away from a per rev value (Af = 0.1 per rev in 
Eq. 6). 

Study on the Effect of Initial Designs 

The IADS multilevel decomposition optimization pro- 
cedure is demonstrated for three examples using the three 
starting points shown in Fig. 8. Example 1 (Fig. 8a) uses a 
rectangular planform with a pretwist of -9.0 degrees, root 
chord of 0.3449 ft, and upper level stiffnesses design vari- 
ables initialized to be consistent with the lower level initial 
wall thickness and lumped areas (i.e., matched stiffnesses). 
Example 2 (Fig. 8b) uses a tapered planform with a 
pretwist of -16.0 degrees, root chord of 0.45 ft, and 
matched stiffnesses. The blade is rectangular to 80 percent 
radius and then tapers linearly to the tip with a 3-to-l taper 
ratio. Example 3 (Fig. 8c) uses the same planform and 
pretwist as Example 2 but the upper and lower level stiff- 
nesses are not matched. All these examples use a value of 
-0.4 for the coordination parameter e in Eq. 26. The impor- 
tance of the choice of e is examined in a later section of the 
paper. 
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Example 1 Rectangular Planform ("Initially 
matched stiffnesses") - The starting point for the opti- 
mization is the rectangular blade shown in Fig. 8a. The 
upper and lower level stiffnesses are matched since the up- 
per level stiffnesses are started with the stiffnesses deter- 
mined by the initial lower level design variables. This is 
an infeasible starting point because the lower level stress 
constraints at the root are violated. Results are given in 
Table 3. The initial and final values for the blade planform, 
performance measures, and dynamics measures are given 
in the upper portion of Table 3. The initial and final values 
for the constrained frequencies are given in the middle por- 
tion of the table. Notice that the final value for the fourth 
bending frequency f j, 4 is in a different frequency range 

than the initial value. Final values for the lower level de- 
sign variables and the upper level stiffnesses are given in 
the bottom portion of the table. The final design is able to 
improve the performance characteristics from the ini tial 
blade and satisfy all the constraints. Compared to the ini- 
tial values, the final design represents a 2.1, 2.3, 2.3, 47.6, 
and 3.2 percent reduction in the power required in hover, 
forward flight, and maneuver; hub shear; and upper level 
objective function, respectively. 

The final stiffness distributions for the upper (required 
values) and lower levels (attainable values) are shown in 
Fig. 9. The matching of the chordwise bending stiffness 
EI XX (Fig. 9a), the flapwise bending stiffness EI ZZ (Fig. 

9b), and the torsional stiffness GJ (Fig. 9c) are extremely 
good. As shown in Fig. 9d, the lower level is able to ob- 
tain an extensional stiffness distribution higher than the 
minimum requirement set by the upper level. 

Convergence histories of the individual terms of the 
lower level objective function (Eq. 13) are shown in Fig. 
10 for the three locations - the root (Fig. 10a), point of ta- 
per initiation (Fig. 10b), and the tip (Fig. 10c). Each term 
(denoted stiffness deviation) is a measure of how well the 
upper and lower stiffnesses match. Initially, the stiffnesses 
are matched, but the stress constraints are violated at the 
root. Therefore, the lower level design variables must 
change to satisfy these constraints while keeping the upper 
and lower level stiffnesses matched as close as possible. 
Notice that the chordwise stiffness at the root, torsional 
stiffness at the point of taper initiation, and flapwise stiff- 
ness at the tip are the last stiffnesses to match. Further, it 
appears that stiffnesses at the point of taper initiation are 
particularly difficult to match. This difficulty may be due 
to the fact that the point of taper initiation is a design vari- 
able while the root and tip positions are fixed. 

The reason for the deviations in the stiffness is that the 
upper and lower levels are in conflict. One component of 
the upper level objective function is the hub shear which 
can be reduced significantly by increasing the blade stiff- 
nesses. On the upper level if the optimizer did not have to 
be concerned with stiffness matching, it would increase the 


upper level stiffnesses. Without the lower level to keep the 
stiffnesses in check, a heavy or nonbuildable blade might 
result. 

The information shown in Fig. 10 is collected and used 
to determine when an upper level design variable move 
limit adjustment is necessary during the overall optimiza- 
tion process (recall that approximate analysis is used on the 
upper level and exact analysis is used on the lower level). 
At the present time no automatic move limit adjustment in 
the approximate analysis on the upper level is used. 
Instead, the IADS procedure is run for eight cycles and 
then the stiffness deviations are examined. When the stiff- 
ness deviation increases (e.g. Cycle 16), the design variable 
move limits are manually reduced and the optimization 
process continued for another 8 cycles. In practical 
applications, the optimization procedure would terminate 
after about 30 cycles. However, for demonstration pur- 
poses the convergence criterion is set to a very small value. 
Both the upper and lower levels have the same tight con- 
vergence criterion on each cycle. Overall convergence of 
the IADS procedure might improve if the convergence cri- 
terion is relaxed initially and then tightened as the opti- 
mization proceeds. 

Example 2 Tapered Planform ("initially matched 

stiffnesses") - The starting point for the optimization is 
the tapered Made shown in Fig. 8b. Initially, the upper and 
lower level stiffnesses are matched since the upper level 
stiffness are determined by the lower level design vari- 
ables. However, this is an infeasible starting point since a 
thin wall theory constraint is violated on the lower level. 
The initial and final values for the blade planform, perfor- 
mance measures, and dynamics measures are given in 
Table 4. The final design is able to improve the perfor- 
mance characteristics from the initial blade. However, the 
hub shear increases from the initial value. 

Fig. 1 1 shows the final stiffness distributions for the 
upper (required values) and lower levels (attainable values) 
for the chordwise bending stiffness (Fig. 11a), flapwise 
bending stiffness (Fig. lib), and the torsional stiffness 
(Fig. 11c). As shown in the Fig., the stiffness matching is 
good, although not a good as in Example 1. The lower 
level is able to obtain an extensional stiffness distribution 
(Fig. 1 Id) higher than the minimum requirement. 

Fig. 12 shows the stiffness deviations versus cycle for 
the three matching locations - the root (Fig. 12a), point of 
taper initiation (Fig. 12b), and the tip (Fig. 12c). Early in 
the optimization process, the flapwise and torsional stiff- 
ness are both unmatched. After Cycle 10, the matchings 
improve and after 25 cycles, all three matchings are good. 
At the tip (Fig. 12c) matching proves to be quite difficult. 
The torsional stiffness is the last to match. The reason for 
this is that the blade initial design is tapered and it is diffi- 
cult to place a thin wall section in the space near and at the 
tip and still match the stiffness required on the upper level. 


10 



Example 3 Tapered planform ("initially un- 
matched stiffnesses") - In the previous examples, the 
starting points used matched stiffnesses. The purpose of 
this example is to demonstrate how the IADS procedure 
behaves when it is started from an inconsistent set of stiff* 
nesses (i.e., unmatched stiffnesses). The starting point for 
the optimization is shown in Fig. 8c. The initial stiffnesses 
used in the tipper level are much larger than the stiffnesses 
obtained from the lower level design variables. The initial 
and final values for the blade planform, performance mea- 
sures, and dynamics measures are given in Table 5. The 
power required for all three flight conditions has increased 
substantially along with the hub shear. The initial and final 
constrained frequencies are also included. Notice that a 
bending frequency fb 6 has shifted frequency intervals. 

The final upper and lower level stiffnesses are shown in 
Fig. 13. As shown in Fig. 13, the optimization procedure 
is able to match the upper and lower level stiffnesses suc- 
cessfully. Fig. 14 shows the stiffness deviations for the 
three matching locations - the root (Fig. 14a), the point of 
taper initiation (Fig. 14b), and the tip (Fig. 14c). As shown 
in the figure, after 25 cycles the optimization procedure is 
able to match all three stiffnesses, but it is at the expense of 
upper level performance (see Table 5). From these results 
it appears that while the optimization procedure will work 
when starting from an initial point which has unmatched 
stiffnesses, it is better to start with a set of consistent stiff- 
nesses. 

Observations on the Effect of Initial Design Study 

The IADS procedure has been exercised for three start- 
ing blade planforms - a rectangular planform with matched 
stiffnesses, a tapered planform with matched stiffnesses, 
and a tapered blade with unmatched stiffnesses. In all 
cases the procedure is able to find converged feasible de- 
signs. Comparing Examples 1 and 2 (Tables 3 and 4, re- 
spectively), the reader will find two different final blade 
designs (i.e., design variable values are different) with es- 
sentially the same objective function value. Apparently, 
there are many different combinations of design variables 
which satisfy the matching conditions and more than one 
local minimum. The final solution depends on initial con- 
ditions. In Example 3 (Table .5), it appears that the opti- 
mizer converges to a suboptimal solution when compared 
with Example 2. Both examples started from the same 
planform, but Example 2 starts with matched stiffnesses 
and Example 3 starts with unmatched stiffnesses. Since 
the initial matching of the stiffnesses is relatively easy, this 
suggests that the initial matching should always be en- 
forced. 

Comparing all three examples, the reader will also no- 
tice that each initial blade has a different frequency range 
for the bending and torsional frequencies and each final 
blade design has a frequency which has shifted a frequency 


interval (e.g., f|, 6 in Example 3). During the approximate 

analysis, the optimizer can change the upper level design 
variables such that a frequency can shift intervals. 
However, as the design variable move limits are reduced, 
this shifting is less likely to occur. 

At the present time no automatic move limit adjustment 
in the upper level approximate analysis is used. However, 
the stiffness deviation information (e.g., Fig. 10) can be 
collected and used to determine when an upper level design 
variable move limit adjustment is necessary during the 
overall optimization process. 


Study on the Effect of the Coordination Parameter c 

The purpose of this study is to demonstrate the effect of 
E in the coordination constraint (Eq. 26) on the optimiza- 
tion procedure. Results for three e values (+0.4, -0.2, and 
-0.4) are presented in Table 6 and Figs. 15-16, and 9, 
respectively. If e is a large positive value, the levels are 
essentially independent. The upper level is free to change 
the upper level stiffness and chord distributions in any way 
which will reduce the upper level objective function. The 
only requirement is that the overall stiffness matching 
should not degrade by more than the amount e from the 
best match found on the last lower level optimization. For 
example if e is 0.4, the stiffness matching can degrade by 
40 percent and still satisfy the coordination constraints. It 
is therefore possible that the procedure could converge 
with the upper and lower level stiffnesses being mis- 
matched by as much as 40 percent. A negative value for e 
means that the upper level must improve the matching 
achieved on the lower level by that amount. This section 
of the paper presents results for several values of e using 
the starting point in Fig. 8a which is also used in Example 
1 . 

One choice for e would be zero. This would mean that 
the upper level cannot degrade the matching achieved on 
the lower level. This value was found to be too restrictive 
for the optimization process and the procedure converged 
in three cycles with very little change in the upper level de- 
sign variables. The reason for this can be seen by examin- 
ing the coordination constraint (Eq. 26). At the start of the 
upper level optimization, the coordination constraint at 
each matching location is active (i.e., g=0) since F u is 
equal to F^ . As the upper level optimizer tries to change 

the upper level design variables, the coordination con- 
straints become violated. Therefore, the upper level opti- 
mizer makes only small changes and the process converges 
in three cycles. 

As shown in Table 6, when e is 0.4, the optimization 
process is able to improve the performance and dynamics 
measures over the initial blade values and improve the 
lower level (satisfy the stress constraints). This improve - 
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ment is achieved at the expense of stiffness matching. Fig. 
15 shows the final stiffness distributions for the upper and 
lower levels. The lower level is not able to find a set of 
stiffnesses to match those required by the upper level. This 
final result is technically a feasible design since all the 
constraints are satisfied. Recall the upper and lower stiff- 
nesses need only be as close as possible (lower level objec- 
tive function). The upper level coordination constraints do 
not require the upper and lower level stiffnesses to match 
exactly. 

When e is -0.2, the optimization procedure is able to 
obtain a design that has some improvement over the initial 
starting point (Table 6). The upper level objective function 
is reduced slightly, but not as much as when t is positive. 
As shown in Fig. 16, the upper and lower stiffnesses match 
well for the chordwise stiffness (Fig. 16a), the flapwise 
stiffness (Fig. 16b), and the torsional stiffness (Fig. 16c). 
The lower level is able to obtain an extensional stiffness 
which is slightly larger than that required by the upper 
level. 

Of the values used in this work, the best value for e is 
-0.4. With this value of e, the optimization procedure is 
able to obtain improvement on the upper level and find a 
set of consistent stiffnesses on the lower level. These re- 
sults (Example 1) are included in Table 6 for completeness. 
The stiffness distributions are shown in Fig. 9. 

As shown above, positive values of e result in upper 
level improvement but poor stiffness matching and nega- 
tive values of e result in both upper level improvement 
(although not quite as good as when e is positive) and good 
stiffness matching. This suggests that a gradual reduction 
from a positive to a negative value for e could be benefi- 
cial. The IADS procedure was run with a value of +0.4 for 
e for 8 cycles, +0.2 for 8 cycles, -0.2 for 8 cycles, and fi- 
nally -0.4 for 8 cycles. This technique of gradually reduc- 
ing the value of e did not work. It is felt that the upper 
level planform area and upper level stiffnesses increased to 
improve the upper level objective function when e was 
positive so that by the time e was negative the stiffness 
matching was achieved at the expense of performance and 
dynamic improvement on the upper level. This situation is 
analogous to Example 3 where the mismatched initial con- 
ditions resulted in stiffness matching at the expense of up- 
per level improvement. 


Concluding Remarks 

An integrated aerodynamic/dynamic/structural (IADS) 
optimization procedure for helicopter rotor blades has been 
developed. The procedure combines performance, dynam- 
ics, and structural analyses with a general purpose opti- 
mizer using multilevel decomposition techniques. At the 
upper level, the blade structure and response are repre- 
sented in terms of global quantities (stiffnesses, mass, and 


average strains). At the lower level, the blade structure and 
response are represented in terms of local quantities 
(detailed dimensions and stresses). 

The IADS procedure consists of an upper level opti- 
mization, a lower level optimization, and a coordination 
task. The upper level objective function is a linear combi- 
nation of performance and dynamic measures. Upper level 
design variables include pretwist, point of taper initiation, 
taper ratio, root chord, blade stiffnesses, tuning masses, 
and tuning mass locations. Upper level constraints consist 
of limits on power required in hover, forward flight, and 
maneuver; airfoil drag; minimum tip chord; trim; blade 
natural frequencies; autorotational inertia; blade weight; 
and average strains. 

The lower level optimization sizes the internal blade 
structure to provide the stiffnesses required by the upper 
level and assure the structural integrity of the blade. The 
lower level design variables are the box beam wall thick- 
nesses and several lumped areas which are analogous to 
longitudinal stringers in a wing box cross section. The 
lower level objective function is a measure of the differ- 
ence between the upper level stiffnesses and the stiffnesses 
computed from the wall thicknesses and lumped areas. 
The lower level constraints are on Von Mises stresses, ex- 
tensional stiffnesses, thin wall theory, and dimensional 
limits. 

The coordination task consists of a set of upper level 
constraints which link the levels and promote consistency 
between the upper and lower level stiffnesses. A coordi- 
nation parameter is included in each constraint. This pa- 
rameter specifies how much the upper level can degrade or 
must improve the overall stiffness matching achieved on 
the lower level and may also be interpreted as a measure of 
how closely-coupled the two levels are. It is found that a 
proper value for the coordination parameter is crucial to the 
success of the IADS procedure. If the parameter has a pos- 
itive value, the procedure will converge but the final stiff- 
ness matching can be unacceptable. If the parameter has 
too small of a value (approximately zero), the optimization 
process will terminate without improving the dynamics or 
performance measures. A small negative value for the co- 
ordination parameter encourages the upper level to im- 
prove dynamics and performance using stiffness values 
which the lower level can match. 

The IADS procedure is demonstrated using a model - 
size rotor blade for several initial blade plan forms and 
varying amounts of coupling between the levels. In all 
cases, the IADS procedure achieves successful results. It 
converges to a feasible design regardless of whether or not 
the initial design had a set of consistent stiffnesses. 
However, initializing the upper level stiffnesses with the 
stiffnesses calculated from the lower level design variables, 
greatly improves the final design. 
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The IADS procedure exploits the couplings and inter- 
actions between the disciplines of aerodynamics, dynamics 
and structures. It provides an efficient method to integrate 
structures and/or structural properties into an optimization 
procedure since it guarantees that a structure with a consis- 
tent set of structural properties can be found. The IADS 
procedure provides an optimization technique that is com- 
patible with industrial design practice in which the aerody- 
namic and dynamic design is performed at a global level 
and the structural design is carried out at a detailed level 
with considerable dialogue and compromise among the 
groups. 
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Appendix A - Coordination Constraint 

In a multilevel decomposition approach, the coupling 
between levels is done through a coordination procedure 
(e.g. Refs. 20 and 25). In the present work, the coordina- 
tion procedure based on Ref. 20 is used to reconcile the 
stiffnesses required on the upper level with the stiffnesses 
the lower level can actually obtain. This reconciliation re- 
sults in one upper level constraint at each matching loca- 
tion 


8 = F U - (1 + e)Fq < 0 


(Al) 


where fJj" is the most recent value of the lower level ob- 
jective function (i.e., optimum value of Eq. 13), F u is an 
estimate of the change in F^ which would be caused by a 

change in the upper level design variable values, and e is 
denoted the coordination parameter. This coordination pa- 
rameter specifies how much the upper level can degrade or 
must improve the overall stiffness matching achieved on 
the lower level and may also be interpreted as a measure of 
how closely-coupled the two levels are. If e has a positive 
value, the two levels are not closely-coupled (i.e., they are 
essentially independent). The upper level can change the 
upper level stiffness and chord distributions in any way 
which will improve the upper level objective function as 
long as the stiffness matching is not degraded by more than 
the amount e. If e has a negative value, the two levels are 
closely-coupled and the upper, level is commanded to im- 
prove the matching by the amount e. 

Eq. Al is the general form of the coordination con- 
straint as formulated in Ref. 15. The form of the coordina- 
tion constraint used in this work is obtained by approximat- 
ing F u in terms of the current optimum lower level objec- 
tive function F^\ If F^ is expanded in terms of a first or- 
der Taylor series about the lower level optimum, then F u 
can be approximated by 


where 


9F l 

9DVi 


is the derivative of the optimum lower level 


objective function with respect to the upper level design 
3KS 

variables, is the derivative of the active lower level 

9DVi 

constraint (Eq. 25) with respect to the upper level design 
variables, X T is the Lagrange multiplier given by 



(A4) 


where — — is the derivative of the active lower level con- 
dv 

straints with respect to the lower level design variables at 

T 

the lower level optimum. At a lower level optimum, X 

T 

will be positive. If no lower level constraint is active, X 
is set to zero. Substituting Eq. A2 into Eq. Al, the coordi- 
nation constraint g is approximated by 




f‘- + T dpL 

ADVj 

° it, dDVi 

0 


-[l + e]F^0 (A5) 


or simplifying 


g = 


N D v dF L 


i=l 


dDV 


ADV 


to 


-eFL<0 


(A6) 


Substituting Eq. A3 into Eq. A6, the coordination con- 
straint becomes 


g = 


NDV 

I 

i=l 


3F l 


3D Vi 


ADV_ X T 9K1 

, ^lo 


- eFq <, 0 (A7) 


which is the form implemented in this work. 


13 



The derivative of the coordination constraint is obtained 
by differentiating Eq. A7 with respect to upper level design 
variables 


h 

dDVj 



(A8) 


Appendix B - Lower Level Structural Analysis 

The purpose of this appendix is to summarize the ele- 
mentary equations describing the geometry and structural 
analysis for the lower level structure. A typical cross sec- 
tion of the thin-walled isotropic box section is shown in 
Figure 5. For simplicity the top and bottom wall thick- 
nesses, tj and tj are equal. The total cross-sectional area, 

A is the sum of the cross-sectional areas of the box beam 
elements Aj and the lumped areas aj (described in the main 
text) 


n m 

A=lAi+Xaj (Bl) 

i=l j=l 

Using the familiar relations, the centroid of the cross-sec- 
tion is calculated from the following equations 


X c = 


n m 

s A i x i + i>jXj 

i=l j=l 


(B2) 





12 

k=l n+m 

(B4) 

h k b^ 



12 

k=l,...,n+m 

(B5) 


where is the base of the kth rectangular element and h k 
is the height relative to the X axis and I^. is equal to zero 

for symmetric elements. Using the parallel axis theorem, 
the moments of inertia of each element are found with re- 
spect to the centroid of the box beam as follows 

f\ 

^cxir = + A k4 k 

2 (B6) 

IcZfc - Izfc + A k C k 

where I CXk , and I CZk are the moments of inertia of the 

kth element about the centroid of the box beam, 
I Xfc and I Zk are the moments of inertia of the kth element 

about its centroidal axes, d k and c k are the distances from 

the centroid of the element to the centroid of the box beam 
in the X and Z directions respectively. The total moments 
of inertia for the box beam are equal to the sum of the ele- 
ment inertias. 


I XX “ S^CXfc 

Izz ~ Xlcz k 


(B7) 


The polar moment of inertia for the box beam is calculated 
using the method described in Ref. 38. 


and 


Z c - 


n m 

X A i Z i + 5>jZj 

i=l j=l 


(B3) 


where Xj and Z x are coordinates in the chordwise and 

flapwise directions respectively that specify the distance of 
the centroid of the ith element area A i from the reference X 

and Z axes shown in Figure 5. Similarly, xj and Zj are co- 
ordinates that specify the distance of the centroid of the jth 
lumped area aj from the reference axes, n is the number of 

elements that the cross section is divided into for ease of 
calculations and m is the number of lumped areas. 


Next, the area moments of inertia of each element about 
its centroidal X and Z axes are calculated from 



(B8) 


where A c is the enclosed area of the mean periphery of the 

box beam wall, ds is the differential circumferencial length 
along the box beam, and t is the local thickness of the wall. 

In order to calculate the lower level objective function, 
the bending and torsional stiffnesses of the box beam are 
necessary. For an isotropic beam the moments of inertia, 
I xx and I zz , calculated above are multiplied by Young's 

Modulus E to acquire the bending stiffnesses, EI XX and 
EI^, in the chordwise and flapwise directions respectively. 

Similarly, the polar moment of inertia is multiplied by the 
torsional modulus of elasticity G to acquire the torsional 
stiffness of the beam, GJ. 
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The stresses for the constraints in the lower level opti- 
mization are evaluated at the comers of the box beam us- 
ing Von Mises stress measure which is given by 



V(o,t) = ^o 2 + 3t 2 


(B9) 


the 39th Annual Forum of the AHS. May 9-11, 1983, 
St. Louis, Missouri.. 


where o is the axial bending stress at the outer fiber of the 
cross section 


'=fe)w 


M 


xx 


V ^xx 


CF 

Zouter+— (BIO) 


and T is the shear stress due to torsion in the wall of the 
section with thickness t 


8. Davis M. W. and Weller, W. H.: Application of 
Design Optimization Techniques to Rotor Dynamics 
Problems. Journal of the American Helicopter 
Society, Vol 33, No. 3, July 1988. 

9. Celi, R. and Friedmann, P. P.: Efficient Structural 
Optimization of Rotor Blades with Straight and Swept 
Tips. Proc. of the 13th European Rotorcraft Forum, 
Arles, France, September 1987. Paper No. 3-1. 


T _ m t 

2Agt 


(Bll) 


where is the flapwise moment; M xx is the lag mo- 
ment; CF is the centrifugal force; and M T is the torque at 

the section. The shear stress due to transverse loads has 
been neglected for simplicity. M a , M xx , CF, and M T are 

computed in the upper level analysis for forward flight and 
maneuver, multiplied by a factor of safety and then 

passed to the lower level. 
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Tabl^^^esigr^ariabl^ounds 



Design 
Lower bound 

variables 
Upper bound 

Twist (deg) 

-20.0 

-5.0 

Point of taper initiation (r/R) 0.26 

0.985 

Taper ratio 

0.05 

5.0 

Root chord (ft) 

0.05 

0.833 

EI XX (Ibm-ft 2 ) 

50.00 

20000000.0 

EI^ (lbm-ft 2 ) 

5.00 

1000.0 

GJ (lbm-ft 2 ) 

5.00 

1000.0 

EA (lbm) 

1000.00 

200000000.0 

m, (slug/ft) 

0.0 

0.50 

Xi (r/R) 

0.24 

0.95 

ti (ft) 

0.00008 

0.01 

a i (ft 2 ) 

0.0 

0.00004 


Tabl^^ftrametei^an^lightcon^Uons 

Parameters 


Minimum autorotational inertia. AI 

min 

23.69 lbm-ft 2 

Allowable drag coefficient, 

0.12 

Minimum tip chord, c, 

iimn 


0.083 ft 

Number of blades, N 


4 

Number of aerodynamic segments 


HOVT 


19 

CAMRAD/JA 


18 

Number of structural segments 

50 

Number of design variables 



Upper level 


22 

Lower level 


18 (6 per 
location) 

Power available, P 0 
« 


20 hp 

Blade radius, R 


4.68 ft 

Maximum blade mass, W 


3.5 lbm 

Factor of safety, L 


2.0 

Af 


0.1 per rev 

Allowable average strain, e a 

i 

0.05 ft/ft 

ITER_ 

max 


40 

Pmax 


300 

Allowable stress, a a 


8.352xl0 6 h/ft 2 

Young’s modulus, E 


15.26xl0 8 lb/ft 2 

Flight 

conditions 


Rotational velocity 
Hover tip Mach 

639.5 RPM (in Freon 
density of 0.006 slug/ft 3 ) 

0.628 


number 

Hover 

Forward flight Maneuver 

C L 0.00810 

0.00810 

0.00985 

C D 

-0.000811 

-0.000596 

Advance ratio - 

0.35 

0.30 


Table 3 - Rectangular planform starting point with matched 
stiffnesses - Example 1 




Initial 

Final 

Hover power (hp) 

14.81 

14.50 

Forward flight power (hp) 

13.26 

12.96 

Maneuver power (hp) 

12.22 

11.94 

Hub shear (lb) 

2.1 

1.1 

Objective function 

20.58 

19.91 

Twist (deg) 


-9.0 

-11.47 

Taper initiation 

0.7 

0.701 

Taper ratio 


1.0 

1.66 

Root cbord(ft) 


0.3449 

0.3770 

mj (slug/ft) 


0.0 

0.0002761 

m 2 (slug/ft) 


0.0 

0.003199 

m 3 (slug/ft) 


0.0 

0.002014 

y» (f0 


- 

0.450 

y2 ( ft ) 


- 

0.583 

y 3 (ft) 


- 

0.453 

Cycles to converge 

- 

76 

f b3 per rev 



2.68 

f b4 per rev 


3.77 

4.57 

f b5 per rev 


4.52 

4.88 

fb 6 Per rev 


7.22 

7.55 

f tl per rev 


7.30 

7.30 

f t2 per rev 


3.61 

3.83 

Final Design 

Root 

Point of taper 

Tip 

Variables 


initiation 


Lower Level 




t, (ft) 

0.002366 

0.002427 

0.0004517 

t 2 (ft) 

0.003261 

0.009954 

0.0003766 

t 4 (ft) 

0.003414 

0.009954 

0.0003766 

a, (ft 2 ) 

0.00003341 

0.00003293 

0.00001610 

32 (ft 2 ) 

0.00001615 

0.00003084 

0.00001192 

a 3 (ft 2 ) 

0.00003281 

0.00003293 

0.00001610 

Upper Level 




EI„(lb-f t 2) 

2057.0 

2974.1 

153.7 

EIja (lb-ft 2 ) 

122.21 

140.03 

8.61 

GJ (lb-ft 2 ) 

127.93 

128.53 

5.87 

EA (lb) 

797370. 


212230. 
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Table 4 - Tapered planform starting point with matched 
stiffnesses - Example 2 


Hover power (hp) 

Forward Flight power (hp) 
Maneuver power (hp) 

Hub shear (lb) 

Objective function 

Twist (deg) 

Point of taper initiation 
Taper ratio 
Root chord (ft) 
nt| (slug/ft) 
m 2 (slug/ft) 
m 3 (slug/ft) 

y, (ft) 
y2 (ft) 

y 3 (ft) 

Cycles to converge 

f bj per rev 
f„ 4 per rev 
f b5 Per rev 
f„ 6 per rev 
f t , per rev 
f, 2 per rev 


__Final_ 

14.74 

13.02 

11.84 

0.66 

19.93 

-10.85 

0.37 

1.64 

0.4932 

0.008961 

0.01354 

0.0246 

0.24 

0.616 

0.622 


Table 5 - Tapered planform starting point with unmatched 
stiffnesses - Example 3 


Hover power (hp) 

Forward flight power (hp) 
Maneuver power (hp) 

Hub shear (lb) 

Objective Function 
Twist (deg) 

Point of taper initiation 
Taper ratio 
Root chord (ft) 
mj (slug/ft) 
m 2 (slug/ft) 
m 3 (slug/ft) 

yi (ft) 

Y2 (ft) 

y 3 (ft) 

Cycles to converge 
f bj per rev 
f b+ per rev 
f bj per rev 
fh< per rev 
f tj per rev 
f, 2 per rev 


Final 

16.64 

17.46 

14.89 

2.45 

24,62 

-11.98 

.889 

1.3148 

0.7364 

0.008546 

0.007797 

0.009030 

0.323 

0.439 

0.393 

92 

2.90 
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Table 6- Effect of e on multilevel decompostion 




izatiom 


ure 


Initial Final Rnal_ 

§ +.4 -1 

Hover power 14.81 14.44 14.60 

(bp) 

Forward flight 13.26 12.77 13.11 

power (hp) 

Maneuver 12.22 11.75 11.96 



power (hp) 


Hub shear (lb) 

2.1 

0.2072 

1.85 

1.1 

Objective 

20.58 

19.48 

20.22 

19.91 

function 





Twist (deg) 

-9.0 

-13.32 

-11.12 

-11.47 

Point of taper 

0.7 

0.786 

0.825 

.701 

initiation 





Taper ratio 

1.0 

3.16 

1.41 

1.66 

Root chord (ft) 

0.3449 

0.3651 

0.3606 

0.3770 

mj (slug/ft) 

0.0 

0.02571 

0.00135 

0.0002761 

m 2 (slug/ft) 

0.0 

0.00211 

0.0000995 

0.003199 

m 3 (slug/ft) 

0.0 

0.00099 

0.0000727 

0.002014 

yi (ft) 

- 

0.4124 

0.3115 

0.450 

y 2 (ft) 

- 

0.4154 

0.3950 

0.583 

y 3 (ft) 

- 

0.4382 

0.4292 

0.453 

Cycles to 

- 

90 

152 

76 


converge 



Figure 1. General multilevel optimization procedure with 
two levels. 
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Figure 2. Multilevel strategy for rotor blade aerodynamic, 
dynamic and structural optimization 


L t -1 




lina of 
•ymmatry 

M 

h 

1 

*1 | *2 | *3 

"VI | *2 | *3 

* 

— i 

— *2 

a b ► 

T* 


Figure 5. Lower level design variables 



Figure 3. Upper level analysis flowchart 



Loads, stiffnass, and local chord 
from uppar lava I 


Calculata stlffnassas from 
lowar laval daslgn varlablas 


^ 




— 

Structural Analysis 
Load casa 1 
Forward Flight 


Structural Analysis 
Load casa 2 
Manauvar 












Objactlva function 
and constraints 







| Optimlzar 


jupdatad daslgn varlablas 



no 


Figure 6. Lower level flowchart 


Figure 4. Upper level design variables 
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Figure 7. Overall optimization flowchart 
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Figure 8. Initial designs used to start IADS procedure 
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Rgure 9. Final upper and lower level stiffness 
distributions - Example 1 
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10b) Point of taper initiation 
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Figure 10. Convergence history of upper and lower level 
stiffness deviations - Example 1 . 
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lid) Extensional stiffness 


Figure 11. Final upper and lower level stiffness 
distributions - Example 2 
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12b) Point of taper initiation 
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13a) Chordwise bending stiffness 
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Figure 12. Convergence history of upper and lower level 
stiffness deviations - Example 2. 


Figure 13. Final upper and lower level stiffness 
distributions - Example 3 
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Figure 14. Convergence history of upper and lower level 
stiffness deviations - Example 3. 






15d) Ex tensional stiffness 


Figure 15. Final upper and lower level stiffness 
distributions for e=+0.4 
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16a) Chord wise bending stiffness 
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Figure 16. Final upper and lower level stiffness 
distributions fore=-0.2 
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